Small time heat kernel asymptotics at the cut locus on 

two-spheres of revolution 



Davide Barilari 

CNRS, CMAP Ecole Poly technique, Paris, Equipe INRIA GECO Saclay-ile-de-France, Paris, France 

barilari@cmap .polytechnique . f r 

Jacek Jendrej 
CMAP Ecole Polytechnique, Paris, France 
jacek. jendrej (Spolytechnique . edu 

November 9, 2012 
Abstract 

In this paper we investigate the small time heat kernel asymptotics on the cut locus on 
the class of two-spheres of revolution, which is the simplest class of 2-dimensional Riemannian 
manifolds different from the sphere with non trivial cut-conjugate locus. We determine the 
degeneracy of the exponential map near a cut-conjugate point and present the consequences 
of this result to the small time heat kernel asymptotics at this point. These results give a 
first example where the minimal degeneration of the asymptotic expansion at the cut locus is 
attained. 

1 Introduction 

The study of the heat kernel and the geometric analysis of the Laplace-Beltrami operator on 
Riemannian manifolds is a very old problem and it has been an object of an increasing attention 
during the last century. 

In particular it has undergone a strong development in the last three decades during which 
many results have been obtained relating the properties of the heat kernel (such as Gaussian 
bounds, asymptotics, etc.) to the geometry of the manifold itself (optimality of geodesies, bounds 
on curvature, etc.). A comprehensive introduction to the subject can be found in OUT] and 
references therein. 

In this spirit, the fundamental relation between the metric structure of the manifold and the 
properties of the heat kernel is provided by a result of Varhadan, stating that the main term in the 
small time asymptotic expansion of the heat kernel is given by the Riemannian distance. 

Theorem 1 (see [E]). Let M be a complete Riemannian manifold, d be the Riemannian distance 
and pt its heat kernel. For x,y G M we have 

lim 4t log p t (x,y) = -d 2 (x,y). 

uniformly for (x,y) in every compact subset of M X M. 

Besides this result, which gives the basic estimate for the heat kernel and holds for every pair of 
points on the manifold, one can study how the singularities of the Riemannian distance reflect on 
the small time asymptotics of the heat kernel, to get refined results. By singularities of the distance 
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here we mean the set of points where the Riemannian distance is not smooth. In particular we 
are interested in how the presence of the cut locus (the set of points where geodesies lose global 
optimality) and the conjugate locus (the set of points where the differential of the exponential map 
is not invertible) can affect the asymptotic expansion of the heat kernel. 

This relation has already been established in the case of Riemannian geometry, as stated in 
[2jl 14j,[T6] (see also [TT] for a probabilistic approach to the Riemannian heat kernel) . We recall here 
a corollary of these results that is contained in [2] (where the results are stated in the more general 
framework of sub- Riemannian geometry [UdS]), saying that the leading term in the expansion of 
Pt(x, y) for t — > other than the exponential one depends on the structure of the minimal geodesies 
connecting x and y. 

Theorem 2 (see [2]). Let M be a complete Riemannian manifold, let d denote the Riemannian 
distance and pt the heat kernel. For x,y £ M we have the bounds 

, » c> 



t n/2- - f „-(i/ 2 ) 

for some to > 0, where C\, C*2 > depend on M , x, and y. Moreover 

(i) if x and y are not conjugate along any minimal geodesic joining them, then 

Pt(x,y) = £±^ e -*<*")/« for < t < t . 

(ii) if x and y are conjugate along at least one minimal geodesic connecting them, then 

P^y) > t(TO/2 )+ (1/4 ) e '" 2(X ' y)/4t ^ for0<t< t . (1) 

Note that (i) holds in particular for y close enough to x since the exponential map starting from 
a point is always a local diffeomorphism in a neighborhood of the point itself. 

The goal of this paper is to show some explicit examples of two dimensional manifolds (that 
include ellipsoids of revolution as a particular case) where the inequality (ii) in Corollary [2] is sharp 
at some point y that is both cut and conjugate with respect to x (in the sense that the exponent 
a in the asymptotic pt(x,y) ~ ^—oi e -<f{x,y)/u sa tigfies a = n/2 + 1/4). 

To this purpose, we analyze the small time asymptotic on some particular class of two dimen- 
sional surfaces which are the simplest perturbations of the two-dimensional sphere. Indeed the first 
non trivial example of a complete Riemannian manifold where an asymptotic expansion at the cut 
locus is computed is given by the standard two-dimensional sphere S 2 . In this case the cut locus of 
a point x coincides with the conjugate locus and collapses to the antipodal point x to x. The heat 
kernel on the two sphere was first computed in [7] and an elementary computation shows that the 
small time asymptotics at these points is given by 

p t (x, x) ~ -l^e-^)/ 4 *, for t -> 0. (2) 

Since, at the authors' best knowledge, there is no explicit proof of this expansion, we wrote it 
in Appendix [Aj Notice that in this case n = 2, thus the exponent in the expansion pt(x,y) ~ 
j.-a e -d 2 (x,y)/u j g no ^. minimal one (cf. ([I])), as a consequence of the fact that the cut-conjugate 
point is reached by a one parametric family of optimal geodesies. (See also Remark 1191) 
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The next class of examples that is natural to consider is the one of two-dimensional surfaces of 
revolution, in which of course ellipsoids of revolution are the main model. 

The determination of the cut and conjugate loci on a complete surface, even a two-dimensional 
surface of revolution, is a classical but difficult problem in Riemannian geometry. Even on an 
ellipsoid of revolution, the computation is not a standard exercise. In [3], the foreseen conjugate 
and cut loci were given as a conjecture and the first proof appeared only recently in |12j . For other 
results about the conjugate and cut locus on surfaces of revolution (and some generalizations) one 
can see also [5|il3|[l8] . 

On an oblate ellipsoid of revolution the cut locus of a point different from the pole is a subarc 
of the antipodal parallel. On a prolate ellipsoid it is a subarc of the opposite meridian. In the 
first case the Gaussian curvature is monotone increasing from the north pole to the equator and 
decreasing in the second case. 

This result is also a particular case of a more general one contained in [18] about the cut locus 
of a so-called two-sphere of revolution. 

Theorem 3 (see P2]). Given a smooth metric on S 2 of the form dr 2 + m 2 (r)d9 2 , where r is the 
distance along the meridian and 9 the angle of revolution and m is smooth, assume the following: 

(i) m(2a — r) = m(r) (i.e. reflective symmetry with respect to the equator, where 2a is the 
distance between poles). 

(ii) the Gaussian curvature is monotone non- decreasing (resp. non-increasing) along a meridian 
from the north pole to the equator. 

Then the cut locus of a point different from the pole is a simple branch located on the antipodal 
parallel (resp. opposite meridian). 

In this paper we will deal only with the cut locus of a point that belong to the equator of the two 
sphere of revolution. Due to this restriction we can require the reflective symmetry with respect to 
the equator only locally, that means that our result is stated under the following assumptions: 

(Al) m(r) = ip((a — r) 2 ) in some neighborhood of r = a, where ifi is a smooth function. 

(A2) the Gaussian curvature is monotone non-decreasing along a meridian from the pole to the 
equator. 

Notice that in the condition (A2) we specify the Gaussian curvature to be monotone increasing. 
For instance, in the class of ellipsoid of revolution, this means that we are considering only the 
oblate ones. 

To prove our result about the asymptotic of the heat kernel, we also need a nonsingularity 
assumption of the two-sphere, that is the following: 

(A3) K"(a) where K(tq) is the Gaussian curvature (that is constant) along the parallel r = ro- 

Notice that an oblate ellipsoid of revolution is a typical example of a two-sphere of revolution 
satisfying all the three assumptions (Al), (A2) and (A3). (See also Section [3731 ) 

Using these results, we analyzed the structure of the degeneracy of the exponential map on such 
a two-sphere of revolution around a cut-conjugate point on the equator proving that, under these 
assumptions, the following asymptotic expansion holds. 
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Theorem 4. Let M be a complete two-sphere of revolution satisfying assumptions (Al)-(AS) above, 
and let d be its Riemannian distance and pt its heat kernel. Fix x G M along the equator and let y 
be a cut-conjugate point with respect to x. Then we have the asymptotic expansion 

Pt(x,y)~^e- d2 ^y^, fort^O. 

This result is a consequence of the fact that the degeneracy of the exponential map is in relation 
with the asymptotic of the heat kernel via the asymptotic expansion of the so-called hinged energy 
function (see Section 12.11 for the precise definition and [21 Theorem 24] for the aforementioned 
relation) which in the two dimensional case can be explicitly analyzed. 

1.1 Structure of the paper 

In Section [2] we introduce some preliminary material. In particular Section 12.11 contains the re- 
sults about the relation among the heat kernel asymptotics and the degeneracy of the exponential 
map, while Section [2.21 recalls the formal definition and the basic properties of the two-spheres of 
revolution. Sections 13.11 and 13.21 contain the proof of the main result of the paper. In Section [3731 
these results are expressed in terms of the geometry of ellipsoids. Finally in Appendix [Aj given the 
formula for the heat kernel on S 2 , we compute its expansion at the cut locus, while in Appendix iBl 
we prove some technical lemmas that are needed in the proof. 

2 Preliminary material 

In this section we introduce some preliminary material that is needed for the proof of our main 
theorem. In particular we recall the main results about the relation between the heat kernel 
asymptotics and the degeneracy of the exponential map. For more details one can see [2j. A 
description of these results in the Riemannian case can be also found in [14"1[T6] . 

2.1 Hinged energy function and cut-conjugate points 

In what follows, wherever not specified, M is an n-dimensional complete Riemannian manifold and 
d denotes the Riemannian distance. 

Notation. If x G M and v G T X M we denote by "f v {t) the geodesic starting from x with initial 
velocity v. Moreover we denote by exp : TM — > M the exponential map defined by exp(u) = 7^(1). 
If the initial point is fixed, we use the notation exp x = exp \t x m- 

Definition 5. Let x G M. The energy function from x is the function E x : M — > K defined by 

Ex(y) ■= ^d(x,y) 2 . 

The gradient of the energy function at a point y G M, at a differentiability point, is easily 
computed by means of the unique geodesic joining x and y. We have the following proposition 

(see my- 

Proposition 6. The function E x is smooth in the set V x C M of points that are reached from x by 
a unique minimizing and non conjugate geodesic. Moreover, if y = exp x (v) G V x then VE x (y) = 
7„(1) G T y M. 
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Now we define the main object that is involved in the asymptotic of the heat kernel. 
Definition 7. Let x, y G M. The hinged energy function h xy : M — > K is defined by 

h x>y (z) := E x {z) + E y {z), z G M. 

Observe that the function h xy is smooth on V x n V y . In the next proposition we characterize 
its minima. 

Proposition 8. Let x,y E M and let 7 : [0, 2] — > M be a minimal geodesic from 7(0) = x 
to 7(2) = y. Then m.mh x ^ y = ^d(x,y) 2 and this minimum is attained at the point zq := 7(1). 
Moreover every global minimum of h X)V is a midpoint of some minimal geodesic joining x and y. 

Proof. Let z G M, and set a := d(x,zo), b := d(zQ,y), a := d(x,z), f3 := d(z,y). Then a = b = 
^d(x,y), so h x ^ y (zo) = \a? + \b 2 = jd(x,y) 2 . From the triangle inequality we have a + /3 > d(x,y), 

2 

2 J - I 



thus 



Wz) = ^>F±^ >W- (3) 



Suppose now that the equality in (|3|) holds. Then - \^ = f ^-jjpl and a + /? = y). From this 



2 y 

we deduce that z lies on a minimal curve from x to y and a = /3. Then this curve is a geodesic and 
z is its midpoint. □ 

Remark 9. Notice that the midpoint of a minimal geodesic lies in the "good" set V x (~)V y even if 
y G Cut(x). 

The next proposition shows how the degeneracy of exp^, near y is reflected by the behavior of 
h XjV near the midpoint of a geodesic joining them, which is assumed to lie in the "good" region. 

Proposition 10. Fix x,y G M and let a : (— e, e) -^r V x C\V y <Z M be a smooth curve such that q(0) 
is a critical point of h X:V . Assume moreover that 3 k G N such that Vh Xjy (a(s)) = s k u(s), where u 
is a smooth vector field along the curve a with u(0) / 0. 

Let 7„( s ) be the minimal geodesic joining x to a(s) in time 1 and set f3(s) := 7 V ( S ) (2). 
Then (3 : (— e, e) —¥ M is a smooth function, /3(0) = y and 

lim = kdexp (7,(0) (1)) «(0), (4) 

s— >0 s 

where in the last formula we identify T a ^M with T^ (o) ( 1 )(T Cf ( )M). 

Remark 11. Observe that under our assumptions v(s) is a smooth curve v : (— e, e) — > T X M and 
every such curve v(s) gives rise to the corresponding curve a(s) := 7„( s )(l). 

Remark 12. If a(s) is the curve of midpoints of a one parameter family of optimal geodesies joining 
x to y then /3(s) = 0^ since the curve /3(s) is constant. 

Proof. The first part is obvious, since f3 is a composition of smooth functions. Let 7u>(s) ■ 

R M 

be the optimal geodesic from 7 W ( S )(0) = y to 7w s )(l) = a(s). It follows from Proposition [6] that 

7«( s )(l) = -7W(s)(l) + ^K, y {a{s)) = - j w ( s )(l) + s k u(s). (5) 
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We want to use the Taylor expansion of exp near 7u(o)(l) = ~~ TWo)(l) ^ TM. 

To this end we introduce on TM local canonical coordinates (z, v) £ M. 2n in a neighborhood of 
a ( ) ^ TAf. In these coordinates we write 

7«(0)(1) = -7w(o)(l) = (°' c )> 

7i»(*)(l) = (z(s),c + a(s)), 
-7w(a)(l) = (^(s),c + fe(s)), 
tt(s) = (z(s),u(s)), 

so that a(s) — b(s) = s k u(s). Let 

ex Pi (z, c + /i) = (dexp (0jC) O, h)) + p 2 hf + • • • + tp k (z, h) k + 0(\(z, h)\ k+1 ) 

be the Taylor expansion of the i-ih coordinate of exp near %(o) (1) = (0, c) (here is some j-linear 
map (R 2n y M for j = 2, . . . , k). 

Observe that exp^ (z(s), c + b(s)) = yt and expj (z(s),c + a(s)) = Pi(s). Also, for j = 2, . . . , k 
we have 

because \(z(s),a(s)) - {z(s),b(s))\ = \s k (0,«(«))| = 0{s k ) and \(z(s), a(s))\ + \(z(s),b{s))\ = O(s). 
Hence 

AW = ft + (dexp (0jC) (0, a( S ) - &(*))) . + 0( S fc+1 ) = Vl + s fc (dexp (O)C) (0, «(*))) . + 0(s k+1 ), 
so that 

/3( g ) = A:s fc - 1 dexp (0c) (0,u(s)) +0(s fc ) 
= its i: - 1 dexp (OiC) (0,«(0)) +0(s fc ) 
= ^- 1 dexp^ (0)(1) (n(0)) + O( S fc ), 

and the conclusion follows. □ 

Corollary 13. Xei x,y G M and Zei zq £ M be a critical point of h XtV . Assume additionally that 
zq £ V x r)V y . Then y is conjugate to x along some geodesic if and only if zq is a degenerate critical 
point ofh X)y . 

Proof. Suppose first that y is conjugate to x along 7„(o). Then there exists a smooth curve v : 
(— e, e) — > T X M such that $(0) = 0. Let a : (— e, e) — > M be the corresponding curve of midpoints 
(see Remark II ip . Then in ([4]) we must have k>2. Thus Hess h Xj y(a(0)) (v, Q;(0)) = for any v. 

Conversely, if Hess h Xt y(zo)(v, w) = for all v, it suffices to take an arbitrary regular curve 
a : (— e, e) — > V x D V y such that a(0) = zo and ot(0) = u; to obtain /3(0) = 0. □ 

Remark 14. Hess h x ,y(zo) is never degenerate along the direction of a minimal geodesic connecting 
x and y (where zq is its midpoint). Indeed if 7 : [0, 2] — >• M is such a minimal geodesic, zo is its 
midpoint and a(s) := 7(1 + s), then h Xt y(a(s)) = \d(x, y) 2 (l + s 2 ). 



= 0(s fe+1 ), 
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Let x G M and let 7 : R — > M be a geodesic with 7(0) = x. Suppose that 7(2) = y is a conjugate 
point of x along 7 and that 7|[o,2] is minimal (such a point is called a cut-conjugate point). 

It follows from Proposition [8] and Corollary 1131 that zq := 7(1) is a global minimum of h XtV and 
a degenerate critical point. Assume that M has dimension 2. Then, according to Remark 1141 the 
degeneracy is only in one direction, which allows us to use the following result, called Splitting 
Lemma or Refined Morse Lemma (it is a special case of [101 Lemma 1]). 

Lemma 15. Let f : U C 1" — > R be defined in an open neighborhood of and assume that is 
an isolated minimum of f such that /(0) = df(0) = and dimkerd 2 /(0) = 1. Then there exist 
coordinates (x±, . . . ,x n -i,x n ) around such that f(x\, . . . ,x n -i,x n ) = x\ + ■ ■ ■ + x 2 l _ l + g(x n ), 
where g : R — > R + is smooth and g{z) = 0(z ). 

Remark 16. The function g is not unique, but its order of vanishing at z = is - this is the maximal 
order of vanishing at z = of functions s h-> f(a(s)) — f(zo) for a smooth curve s t-t a(s) such that 
a(0) = Zq. 

Applying Lemma[15]to h XtV , we obtain that there exists a smooth local coordinate system (z±, z-i) 
near zq such that h x>y (zi, Z2) = h XjV (zo) + z\ + g(z2), where g : R — >■ R+ is a smooth function and 
g(z) = 0(z 4 ). Assume that g vanishes at z = with finite order k + 1 G N (thus k > 3). 

We define in these coordinates the smooth curve a(s) := (0, s). It is immediate that Vh x>y (a(s)) 
= s k u(s) for some smooth vector field u with u(0) 7^ 0. Conversely, there are no smooth curves a(s) 
such that Vh Xi y(a(s)) = s k+1 u(s) with a smooth vector field u. Combining this with Proposition 
[TUland Remark [TT] we obtain the following result. 

Corollary 17. Assume that M has dimension 2. Let x,y G M and assume that y is a cut-conjugate 
point of x along 7 : [0, 2] — > M. Then there exists a smooth curve v : (— e, e) —> T X M such that 
7»(o) = 7 and 

d( Ms) (2),y)=0(s k ), 

where k + 1 is the order of vanishing at of the function g described above. It is impossible to 
obtain a higher order of vanishing. 

This is the aforementioned relationship between the hinged energy function and the order of 
degeneracy of the exponential map. We will show that for two-spheres of revolution satisfying 
assumptions (A1)-(A3), the order of degeneracy is lowest possible, that is k + 1 = 4. 

We now briefly describe how the order of vanishing of the function g is related to the small-time 
asymptotics of the heat equation on M. 

Let (M, g) be a complete orientable n-dimensional Riemannian manifold. We denote dvol the 
Riemannian volume form associated with g and compatible with the orientation (which means it 
equals 1 on a positively oriented orthonormal frame). 

The Laplace- Beltrami operator on M is defineco as Af := div(V/), / G C°°(M). The heat 
operator e tA : L 2 (M) — > L 2 (M) is defined by the formula 

e tA f( P )= [ p t (x,y)f(y)dvol(y), f G L 2 (M), 

JM 

1 recall that the divergence divX of a vector field X is the unique function satisfying Lx dvol = (divX) dvol, Lx 
being the Lie derivative. 
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where pt(x,y) G C°°(IR + x M x M) is the heat kernel, i.e. the unique smooth function satisfying 
the following conditions. 

(dt-A x )p t (x,y) = 0, (6) 
lim / p t (x,y)f(y)dvol(y) = f(x), V/ € C°°(M). (7) 

Recall that the heat kernel exists and is unique. Moreover it satisfies pt(x,y) = pt(y,x) > 0. 

There are classical results of the Hodge theory stating that if M is compact then the operators 
A and e tA are simultaneously diagonalizable and the proper vectors are smooth functions. Their 
spectra determine the long-time behavior of the heat flow and contain topological information (see 
for instance [31l9lfT7] ) . 

Here we investigate the short-time behavior, which turns out to be connected to the geometry 
of the manifold. 

Theorem 18. Let x,y G M and assume that there exists only one minimal geodesic joining x 
and y. Let zq be the midpoint of this geodesic and assume that there exists a coordinate system 
(z\, . . . , z n ) near x such that 

K,y(z) = h x , y (z ) + Z^ + • • • + Z> + 0(N 2mi + • • • + |Zn| 2m "), 

for some integers 1 < m\ < • • • < m n . Then there exists a constant C > (depending on M , x 
and y ) such that 

C + o(l) / d(x,y) 2 \ 
Pt(X ' y) = ^E^ eXP l"^^J- 

A full proof can be found in [16], whereas the main ideas are already in [14J. The result was 
extended to sub-Riemannian manifolds in [2]. 

Remark 19. As explained in [21 Remark 2], if there exists a one-parameter family of minimal 
geodesies joining x and y, the theorem is still valid, but it should be understood that some nii is 
infinite. For example let M = S 2 be a sphere of radius 1 and let x, x be two opposite poles. Then 
we have mi = l,m2 = oo, so the asymptotics is (see Appendix |A]) 

In particular, let us return to the situation of Corollary 1 171 We obtain the following result. 

Corollary 20. Let M be a 2-dimensional orientable Riemannian manifold. Let x,y 6 M and 
assume that there exists only one minimal geodesic joining x and y. Assume further that there 
exists no smooth curve v : (— e, e) —> T p M such that 7„(o) : [0, 2] — > M is the minimal geodesic from 
x to y and 

d(7v(s){2),y) =o(s 3 ). 
Then there exists a constant C > ( depending on M, x and y ) such that 

C + o(l) ( d(x,y) 2 \ 
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Figure 1: Picture of the cut locus (blue) and conjugate locus (red) from x on a surface. 
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Figure 2: Notational convention for a 2-sphere of revolution. 



2.2 Two-spheres of revolution 

In this section we introduce the main class of two dimensional manifolds that we are going to study 
in the sequel and we recall some basic properties. 

Definition 21. Let (M,g) be a compact Riemannian manifold homeomorphic to a 2-sphere. M is 
called a 2-sphere of revolution if there exists a point p E M called a pole such that for any q\ , q<i G M 
satisfying d(p, qi) = d(p, 92) there exists an isometry / : M — > M satisfying f(qi) = q2, f(p) = P- 

Let (M, g) be a 2-sphere of revolution and let p be a pole. It can be proved that p has a unique 
cut point q which is also a pole [181 Lemma 2.1]. Therefore M \ {p, q} can be parametrized by 
geodesic polar coordinates around p, which we denote (r, 9). Let M' := M \ {p, q}. This allows 
to express the Riemannian metric on M' as ds 2 = dr 2 + m(r) 2 d9 2 where mis a positive function 
satisfying lim^o m(r) = P p. 287, Proposition 3]. Geometrically, m is the distance from the 
rotational axis (see also Figure which shows a section of M by a plane containing the rotational 
axis). 

We denote a := \d{p,q) and b := m(a). The set {(r,9) G M : r = 0} is called the equator. 
Each set {(r, 9) £ M : r = const} is called a parallel. Each set {(r, 9) £ M : 9 = const} is called a 
meridian. 

Fix a point x on the equator and, without loss of generality, assume 8(x) = 0. Next, we recall 
a standard result on geodesies on surfaces of revolution. 

Proposition 22 (see [6J). Let ^(t) = (r(t),9(t)) be a unit speed geodesic on M' . There exists a 
constant v , called the Clairaut constant of ~f, such that 

m{r{t)) 2 8(t) = v. (8) 

Remark 23. Observe that m(r(t))9(t) = cosrj(t), where rj(t) is the angle between j(t) and J$| 7 ( t y 
The constant v can be interpreted as the angular momentum of the geodesic. 
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Proposition 24. Fix y£M. Then we have the following properties: 

1. If \u\ > b, then no geodesic on M' emanating from x satisfies (|8j). 

2. If \v\ = b, then there is exactly one geodesic 7 : M — > M' emanating from x satisfying ([8]). Its 
image is the equator. 

3. If < \u\ < b, then there is exactly one geodesic 7 : K — > M' emanating from x satisfying ([8]) 
and r'(0) < 0. 

For i/£R this unique geodesic will be denoted by c u . 

Proof. We assumed 7(f) to be unit speed, that is r(t) 2 + m(r(t)) 2 6(t) 2 = 1. Thus ([8]) is equivalent 
to 

m(r(t)) z 

which in particular gives 

r(0f = (10) 

1. If \v\ > b, we clearly get a contradiction. 

2. If \v\ = b, we get r(0) = and 9(0) = ^. There is a unique geodesic with initial tangent 
vector (r(0),#(0)) = (0, ^-). By symmetry this geodesic will not quit the equator. By ([8]) 
9 (t) is constant, so the geodesic covers the whole equator. 



3. If \v\ < b, we obtain r(0) = — y fe2 fe / 2 (because we chose r(0) negative) and 6(0) = p. This 

determines the unique geodesic with Clairaut constant v. We have \m(r(t))9(t)\ = \ cosr](t)\ < 
1, so m(r(t)) > \v\ > 0. This ensures that the geodesic does not meet the poles. 

□ 

Assumptions. In what follows we make the following assumptions: 

(Al) m(r) = ip((a — r) 2 ) in some neighborhood of r = a, where Tp is a smooth function. 

(A2) the Gaussian curvature is monotone non-decreasing along a meridian from the pole to the 
equator. 

(A3) K" (a) ^ where K(ro) is the Gaussian curvature (that is constant) along the parallel r = r$. 

The assumption (Al) means that the metric is invariant by reflection with respect to the equator, 
in a neighborhood of the equator itself. 

From (A2) and the Gauss-Bonnet Theorem it follows that the Gaussian curvature on the equator 
is strictly positive. Thus from Lemmas 2.2 and 2.3 in [18] one can obtain that m is strictly increasing 
on (0, a]. For z/, < \v\ < b, we will denote by R = R(v) the unique R 6 (0, a] such that m(R) = 
It is the minimal geodesic distance of c v (t) from the pole p. 
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Theorem 4.1 in [18] states that the cut locus of x is a subset of the equator. For < \v\ < a 
the cut point along c u is the first point of intersection of c u with the equator. This point is given 
by the formula [181 P- 385] (r, 9) = (a, 93(1/)), where 

V(v):=2 ( \ / ( \2 2 - 

J r m{r)^J m[r) z — v A 

From now on we will assume that v > (the case v < is symmetric). It can be shown 
that ^{f) is non-increasing for v £ (0,6) \18\ Lemma 4.2]. Thus the cut-conjugate point of x has 
coordinates (r, 9) = (a,lim u ^ b - ip(v)). We note 9 cut := lim,,^- <p{v) and i cut := b9 cut , so that 
Q,(t C ut) is the cut-conjugate point. 

Remark 25. The geodesic c v starting from x is determined by each of the parameters v, R or rj. 
We recall the relationships between these parameters: v = m(R) and cost/ = jr. In particular b — v 
is of order rf as r] — > 0. 

3 Proof of Theorem |4] 

In what follows, we prove Theorem U] in two steps. First we compute the asymptotic expansion of 
the function <p and then we apply this result to investigate the geodesic variations of the equator. 

3.1 Expansion of (p(is) 

We will now derive the asymptotic expansion of (p(v) for v close to b from the asymptotic expansion 
of m(r) near r = a. Let 

^(s) = b-as + Ps 2 + 0(s 3 ), 
be the Taylor expansion of ip, so that 

m(r) = rp((a - rf) = b - a(a - rf + /3(a - r) 4 + O ((a - rf) . 

Proposition 26. The following asymptotic expansion holds 

^^ + ^w (6 ^ >+0((6 -'' )2) ' <12> 

As an immediate corollary we recover the formula for the cut point along the equator. 
Corollary 27. The coordinates of the cut point along the equator satisfy 9 cnt = J^—^ - 

Remark 28. The Gaussian curvature K of a surface of revolution is constant on parallels, and is 
expressed as K(r) = — (see [SI P- 162]). In particular K{a) = — ^jb^ = ^ on the equator 
(from which one gets also a > 0) and one can recover the first conjugate point along the equator 
via the Jacobi equation. Indeed it is a general fact that the first conjugate point along 7 is 7(t C onj), 
where t CO nj is the first positive zero of the solution of the differential equation 

u(t) + K(j(t))u(t) = 0, u(0) = 0, u(0) = 1. 
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In our case K is constant and the solution of this equation is given by 

FT f [2a\ 

whose first positive zero is i cut = t C oni = m ■ On the equator dt = bd9, so we get indeed 

J v 2oo 

O ^cut 7T 

° cut ~ b ~ 

The proof of Proposition [26] requires two elementary lemmas, whose proofs are postponed to 
Appendix [Bj 

Lemma 29. Let U C M. 2 be a neighborhood of (0, 0) and let f : U — )■ R be a smooth function. Let 
f(x, y) = Yli+j<n-l a ij xl V^ + 0((x 2 + y 2 ) n ^ 2 ) be its Taylor expansion. Then for y > the function 
defined by 

F(0) = 5/(0, 0), F(y) = f{x f dX for y > 0, 

Jo v y - £ 



= & + 6 1 y + --- + 6 n _ 1 y"- 1 + O(y n ), where b k = "Pa^ ^— (13) 

rr; Jo VI -it 2 



is smooth and satisfies 

rl u-? du 

i=o 

Lemma 30. Lei 5 : (— e, e) — )• R fee a smooth function. Let f : {—e,e) 2 — >• R &e defined as 

( g(x) - g{y) ., , 

[ c/'(x) ifx = y. 

Then f is a smooth function. Lf g{x) = ^2k=o a k xk + 0(|£| n+1 ). T/ien 

n 

/(x,y)=J> Yl ^ +0{{x 2 +y 2 rl 2 ). 
k=i i+j=k-\ 

Proof of Proposition \2b\ Recall that v = m(R) and R — > a~ as v — > b~ . Define R := a — R. We 
have 

2m{R) dr 



R 



m(r)^m(r) 2 — m(i?) 2 



2m(a - R) df ^ 



m(a — r)\J m(a — r ) + m(a — R)\J m(a — f ) — m(a — R) 
Being b > it is clear that the map 

2m(a - R) 



(?, R) H- 



m(a — r)\J m(a — r ) + m(a — R) 
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is a smooth function of two variables in a neighborhood of (r, R) = 0. Its Taylor expansion can 
be easily computed from the Taylor expansion of m. As for (m(a — r ) — m(a — i?)) -1 / 2 , it can be 
written as 

1 / R 2 -r 2 



(recall that m(r) = ip((a - r) 2 )). By Lemma[30] we know that — ~l Z} — - is a smooth function of 
(r 2 ,R 2 ) and 



R 2 -f 2 



We have a > 0, so A / — ^ r ^ ON is also a smooth function in a neighborhood of (f, i?) = 0. 

Summing up, the right hand side of ()14[) has the form required in Lemma [29l Performing 
explicitly the computation described above (it is long but straightforward) we get the expression 



<p{y) 

with 



a f(?,R)d? 
o v^R 2 " 



Using the fact that -^==5 = §> an d Jo* ^ d "n = f ' we OD tain the coefficients 

7T 7T 7T 7T (66/3 — « 2 ) 7T 

5 ° = 2 a °° = 62 = 2 a ° 2 + I 020 = 8v/2(6a)3/ 2 " 

In other words we have the expansion 

It follows from the fact that a > that (a — i?) 2 is a smooth function of v = m(R) in a 
neighborhood of v = b. An explicit computation gives 

(a-Rf= b -^ + ^l + 0((b-uf), 

a a A 

which together with (fT5|) leads to (fT2|) . □ 

Remark 31. From this formula it is clear that the case 66/3 = a 2 is going to be singular. It is easy 
to compute that this is equivalent to K"{a) = 0, where K{tq) is the Gaussian curvature on the 
parallel r = tq. 
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3.2 Variations of the optimal geodesic 

Recall that the geodesic c&(i) follows the equator and reaches the cut-conjugate point (r, 9) = 
(a, cu t) for t = t cut . We are now interested in variations of this optimal geodesic. Such a variation 
is given by a smooth curve in the space of parameters, that is a smooth curve s i— > (rj(s),t(s)) £ R 2 
such that 77(0) = and t(0) = t cnt - Recall that n(s) = arccos is the angle between (0) and 

Proposition 32. There exist variations of the optimal geodesic Cb such that 

d(c u(s) {t(s)),(a,9 CVLt )) =0(s 3 ). 

If the two-sphere is nonsingular (that is K"(a) 7^ 0), then there exist no variation of the optimal 
geodesic Cb such that 

d(c u ( s )(t(s)),(a,9 cut j) =o{s 3 ). 

Lemma 33. For < v < b let (r u ,9 cu t) be the first point of intersection of c u with the meridian 
P:={9 = 9 cut }. Then r v = a- ^g#V + Otf). 

Proof. We assume that c„(i) reaches r = R before it reaches 9 = # cu t, which means it goes down 
after intersecting with P (this will be true for 77 small enough). This means that r(i) > for 
9 G [#cut, V 9 ( zy )]- Thus from ([9]) and the Clairaut relation we get 



m = vgrpg^g 9 {t) 



m(r{t)) m(r(t)Y z 
this permits to compute 



dr fit) m(r(t))^m(r(t)) 2 — v 2 
~d9~9{t)~ v : 

u 2 m'(r(t)) 



r\t) 



9(t) 



m(r(t)) 3 ' 



2v^m{r(i)) 2 — v 2 m! (r(t)) 



m(r(t)) 4 

d 2 r _ f\t)9{t) - f{t)9{t) _ (2m(r(t)) 2 - v 2 ) m(r(t))m' '(r(t)) 
~d0 2 ~ ~ (9(t)) 3 ~ v* ' 



In particular 



and 



dr \ b\/b 2 — v 



2 



d9 ) 



d 2 r 



6tan(ry), 



0{R) = 0{n). 
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Thus 



-* = L uO/^I. (uL-/» Ur u 

btan(r,)(<p(iy) - 9 cut ) + / / O(r/)d0d0i. 



Recall that b — v = 6(1 — cos rj) = ^- + 0(?? 4 ). The second term is 0{rf) because we integrate twice 
on intervals of length 0(b — u) = 0{rj 2 ). The conclusion follows by substituting b — v= ^3- + 0(ry 4 ) 
in Proposition [26j □ 

Lemma 34. T/ie distance between (a, C ut) anc ^ the geodesic segment c u ([0; <p( u )]) equals ^~^/^s75~~ 
r/ 3 + 0(r/ 4 ). 



Proof. Let g be a point on the geodesic segment under consideration. 

K, - ,; ,/ 

, — /-J , I — /Of ^ Ronoaf mrf "i"V»o pnm nnf citinn frmn 

with 9q instead of 9 cu t gives 



Suppose that d(q, (a, 9 cut)) = 0(?? 3 )- Meridians are minimal curves, so we obtain a — r-q = 0(t] 3 ). 
Hence, by the triangle inequality, \9q — 9 cut \ = 0(r] 3 ). Repeating the computation from Lemma [331 



(6/36 - a 2 ) Vbir , , , s 

a ~ r ~*- WfaU " +0( " ) - 

Thus 

u (6/36 - q 2 ) V6vr „ , , N 

^K^ ut ))> 1 ; 6 ^J /2 ^ +0 ^) 



□ 



Proof of Proposition \3%\ The first statement is a direct consequence of Lemma [ 

Now let (jj(s),t(s)) be a variation of the equator. By the Gauss Lemma we get t'(0) = 0. Hence 
rj(s) ~ s. The conclusion follows from Lemma [341 □ 

Theorem follows now from Proposition [32] and Corollary 1201 



3.3 Oblate ellipsoid as a two-sphere of revolution 

Let M be an ellipsoid with semi-axes 6, 6, c, where b > c. We denote p, q its northern and southern 
pole respectively. Let a be the distance from the equator to a pole, so that a and 6 have the same 
meaning as before. We will express the expansion of the function r i— > m{r) and v i— > <p{y) in terms 
of the semi-axes 6, c. We will see that an oblate ellipsoid which is not a sphere is nonsingular. 

Proposition 35. Let M be an oblate ellipsoid with axes 6, 6, c, where b > c, is a two-sphere of 
revolution which satisfies assumptions (Al), (A2) and (A3). The function m has an expansion 

m(r) = b-a(a- r) 2 + /3(a - r) 4 + O ((a - r) 6 ) (16) 
where a = and /3 = ^ 4& 24 I c - . 
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Proof. Clearly M is a two-sphere of revolution satisfying assumptions (Al) and (A2). We are left 
to prove Eq. (I16|) and the fact M is nonsingular. 

Let x be a point in the northern half of M and let v be its distance from the rotational axis. 
Let R be the geodesic distance from x to p. Then m(R) = v. On the other hand, a — R is the 
length of an arc of an ellipse, which can be expressed by means of an elliptic integral. In this way 
we obtain 

a - R= L \l-^ + w-fi dt 

Let Z = y/b — v. After substitution t = b — z 2 and some operations on power series the integral 
above transforms into 



f z cV2 46 2 - 3c 2 z 2 

/ ~^ + - — + 0(z 4 )dz, 

Jo Vb 2V2tf/ 2 c V ' 



which results in the expansion 

V2c (4b 2 -3c 2 ) , , - x 
a-R= V 2 + V r- , ' Z z + OX 5 . 

By the Implicit Function Theorem Z is a smooth function of R in a neighborhood of R = a. It 
follows from the expression above that 

Vb~^ =Z= ^-(a -R)- — (4b ' ~ 3C2) (a -Rf + O ((a - Rf) , 



so finally, squaring both sides, 

m{R) = „ = 4 _ + - ^ + o ((„ - Rf). 

This proves f)16[) . Finally, notice that the singularity condition reads (see also Remark 13 ip 

6bf5 — a 2 = = = if and only if b = c. 

c b 

Hence spheres are the only singular oblate ellipsoids of revolution. □ 

A Heat kernel asymptotics on S 2 

Let x G S 2 be the north pole and x £ S 2 be the south pole of a standard sphere of radius 1. After 
suitable renormalizations we obtain from formula (1) in [7] 

p t (x, x) = — ("l) n (2^ + l)e~ {n2+n)t . (17) 



8vr 

n=— oo 



Define f(y) = 8f3 / 27r i/2 e y ^ 4t iye ty l 2 . Using standard rules of computing Fourier transforms it is 
easy to see that 



<f(n) :=-!- / f(y)e- m ydy = ^(2n + l)e- 



(n 2 +n)t 
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By Poisson Summation Formula we know that Y^=-oo <-p{n)e iny = Y^k=-oo fiv + 2&7r). Using this 
identity for y = tt we obtain 

oo oo 

p t (x,x)= £ (-l)>(n)= J2 /((^ + 1 K)- 

n=-oo k=— oo 

Notice that for an asymptotic expansion only two terms in the sum on the right are significant. 
Hence we get 

Pt (x,x) * /(-*) + m = ~ ^i<r*w. 

B Proof of Lemmas 

In this section we give a proof of Lemma [29] and Lemma [301 

Proof of Lemma [HKh The function G(y) = Jq ^Mj^dt is smooth on (— e, e). For y > the change 
of variables x = ty gives -F(y) = G{y). Moreover, by definition F(0) = G(0), hence i 7 is smooth on 

[0,6). 

We have f(ty,y) = c (t) + a(t)y + ... + c n ^{t)y n - 1 + 0(y n ), where c k (t) = T,j=oHk-jV- 
Integrating with respect to t we obtain F(y) = G(y) = bo + b\y + . . . + b n -iy n + 0(y n ), where 
bk = Jo ^=^f ; which is exactly the coefficient bk defined in the statement of the Lemma. □ 

Proof of LemmalM Define h : [0, 1] x (-e,e) 2 -)■ R by x, y) = + (1 - t)y). Then /i(t, •) is 
a continuous family of smooth functions. We have 

f(x,y) = / h(t,x,y)dt, 
Jo 

which proves that / is smooth. The Taylor expansion of g is found by comparing the coefficients 
on both sides of the equality (x — y)g(x, y) = f{x) — f(y). □ 
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